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Abstract —Cell Switch-Off (CSO) is recognized as a promising 
approach to rednce the energy consumption in next-generation 
cellular networks. However, CSO poses serions challenges not 
only from the resource allocation perspective hnt also from 
the implementation point of view. Indeed, CSO represents a 
difficult optimization problem due to its NP-complete nature. 
Moreover, there are a number of Important practical limitations 
in the implementation of CSO schemes, such as the need 
for minimizing the real-time complexity and the number of 
on-off/off-on transitions and CSO-indnced handovers. This arti¬ 
cle introduces a novel approach to CSO based on multiobjective 
optimization that makes nse of the statistical description of the 
service demand (known by operators). In addition, downlink and 
uplink coverage criteria are Included and a comparative analysis 
between different models to characterize Intercell interference 
is also presented to shed light on their impact on CSO. The 
framework distingnishes itself from other proposals in two ways: 
1) The number of on-off/off-on transitions as well as handovers 
are minimized, and 2) the computationally-heavy part of the 
algorithm is executed offline, which makes its implementation 
feasible. The results show that the proposed scheme achieves 
substantial energy savings in small cell deployments where service 
demand is not uniformly distributed, without compromising the 
Qnality-of-Servlce (QoS) or requiring heavy real-time processing. 

Index Terms —Cellular networks, energy efficiency, cell switch- 
off, CSO, multiobjective optimization, Pareto efficiency. 


I. Introduction 

F uture hyper-dense small-cell deployments are expected 
to play a pivotal role in delivering high capacity and 
reliability by bringing the network closer to users 11. How¬ 
ever, in order to make hyper-dense deployments a reality, en¬ 
hancements including effective interference management, self¬ 
organization, and energy efficiency are required El. Given that 
large-scale deployments composed of hundreds or thousands 
of network elements can increase the energy consumption 
substantially, the need for energy efficiency (green communi¬ 
cations) has been recognized by the cellular communications 
industry as an important item in research projects and stan¬ 
dardization activities Em. 

Initial attempts to improve the energy efficiency in cellular 
networks were oriented towards minimizing the power radiated 
through the air interface, which in turn reduces the electro¬ 
magnetic pollution and its potential effects on human health. 
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However, most of the energy consumption (between 50% to 
80%) in the radio access network takes place in base sta¬ 
tions (BSs) in and it is largely independent of the BSs’ load. 
Since cellular networks are dimensioned to meet the service 
demand in the busy hour (i.e., peak demand), it is expected 
that, under non-uniform demand distributions (both in space 
and time), a substantial portion of the resources may end up 
being underutilized, thus incurring in an unnecessary expen¬ 
diture of energy. The problem may become worse in many 
of the scenarios foreseen for 5G, presumably characterized by 
hyper-dense small-cell deployments, hierarchical architectures, 
and highly heterogeneous service demand conditions El. 
Therefore, the idea of switching off lightly loaded base stations 
has been considered recently as a promising method to reduce 
the energy consumption in cellular networks. This framework 
is referred to as Cell Switch-Off (CSO) and it is focused on 
determining the largest set of cells that can be switched off 
without compromising the Quality-of-Service (QoS) provided 
to users. Unfortunately, CSO is difficult to carry out due to 
the fact that it represents a highly challenging (combinatorial) 
optimization problem whose complexity grows exponentially 
with the number of BSs, and hence, finding optimal so¬ 
lutions is not possible in polynomial time. Moreover, the 
implementation of CSO requires coordination among neighbor 
cells and several other practical aspects, such as coverage 
provision and the need for minimizing the number of (induced) 
handovers and on-off/off-on transitions. In practice, optimizing 
the number of transitions, as well as the time required for them, 
is advisable because switching on/off BSs is far from being 
a simple procedure, and indeed, this process must be gradual 
and controlled Ellll. Moreover, a large number of transitions 
could result in a high number of handovers with a potentially 
negative impact on QoS El. 

Although CSO is a relatively young research topic, a 
significant amount of contributions has been made. Hence, an 
exhaustive survey is both out of the scope and not feasible 
herein. Instead, a literature review including, in the opinion of 
the authors, some of the most representative works is provided. 
Thus, in the comparative perspective shown in Table the 
following criteria have been considered; 

• CSO type / architecture: CSO solutions can be classified 
as ‘snapshot’ or ‘traffic profiling’ CSO depending on the 
approach followed to take the on/off decisions. In snapshot 
CSO (e.g., oniihiiniiii), decisions involve the analysis of 
discrete realizations of users, i.e., whenever a CSO decision 
is required, information of every single user in the network 
needs to be available at a central unit where an heuristic or 
optimization procedure is performed. Given its nature, this 
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TABLE I 

Summary of Related Work. 


Ref 

Year 

CSO type / architecture 

Realistic 

ICI model 

Coverage aspects 

Complexity and 
feasibility analysis 

13 

2016 

Snapshot / Semi-distributed 

Full-load 

Downlink 

/ 


10 


2016 

Traffic profiling / Centralized 

Full-load 

Downlink 

Partially 


11 


2016 

Traffic profiling / Semi-distributed 

Load-coupling 

Downlink 

Partially 


12 


2015 

Snapshot / Centralized 

Full-load 

Downlink 

Partially 


13 


2015 

Snapshot / Centralized 

Full-load 

Downlink 

Partially 


14 


2014 

Snapshot / Centralized 

Constant 

Downlink 

X 


15 


2014 

Traffic profiling / Semi-distributed 

Full-load 

Downlink 

/ 


16 


2014 

Snapshot / Centralized 

Constant 

Downlink 

X 


17 


2013 

Snapshot / Centralized 

Full-load 

Downlink 

X 


18 


2013 

Traffic profiling / Centralized 

Load-coupling 

Downlink 

X 


19 


2013 

Snapshot / Centralized 

Constant 

X 

Partially 


20 


2013 

Snapshot / Centralized 

Full-load 

Downlink 

Partially 


21 


2013 

Snapshot / Centralized 

Full-load 

Downlink 

X 


22 


2012 

Snapshot / Centralized 

Constant 

X 

X 


23 


2012 

Snapshot / Centralized 

Full-load 

Downlink 

X 


24 


2011 

Traffic profiling / Semi-distributed 

Full-load 

Downlink 

Partially 


25 


2010 

Snapshot / Centralized 

Full-load 

Downlink 

X 


26 


2010 

Snapshot / Centralized 

Constant 

Downlink 

X 


27 


2010 

Traffic profiling / Semi-distributed 

Full-load 

X 

X 


28 


2009 

Snapshot / Both 

Full-load 

X 

X 


type of solution provides the ultimate performance in terms 
of energy savings. However, feasibility becomes a serious 
issue as it will be shown later on. On the other hand, in 
traffic profiling CSO (e.g., lfT5ll24ll27l l. a certain knowledge 
about the service demand behavior is assumed. 

• Realistic ICI model: Switching off/on BSs modifies not 
only the levels of Intercell Interference (ICI), but also the 
resulting load coupling j^ . The load of each BS can be 
understood as the fraction of resources that are being used 
to satisfy a certain service demand. There are three different 
ICI models that can be used in CSO: 

1) Constant: The assumption is that ICI is always constant 
or it does not exist (e.g., lfT4l [191). In such cases, ICI 
levels are calculated as if all cells were always active no 
matter their actual state (on or off). 

2) Full-load (FL): ICI is only created by active BSs, but 
assuming that they are fully loaded (e.g., 13 da El). 
This approach is reasonable in CSO since traffic is 
concentrated in a subset of BSs which tend to be highly 
loaded. Thus, active BSs always cause ICI to its neigh¬ 
bors while BSs in sleeping mode do not. 

3) Load-coupling (LC): The ICI created by each BS is pro¬ 
portional to its load (e.g., imdsi). This model is more 
realistic and accurate but involves more complexity. 

• Coverage: It indicates whether coverage aspects (in down¬ 
link, uplink, or both) are considered. This is an important 
criterion in CSO to avoid coverage holes. 

• Complexity and feasibility: It indicates whether a complex¬ 
ity and/or feasibility analysis is provided. For instance, in 
case of CSO, how much information is required to be 
exchanged among cells is an important metric. Moreover, 
the number of transitions and handovers in a dynamic 
environment must also be estimated. 

As it can be seen from Table snapshot-based schemes 
require centralized operation, and due to the NP-complete 
nature of the problem, heuristics are the preferred way to 


deal with it. Representative examples of this type of solutions 
include mini Hal i26i where the main idea is, in general, to 
sequentially switch 1) on highly loaded cells (cells that can 
get more traffic), or 2) off lightly loaded ones (cells whose 
users can easily be migrated to neighbor BSs) until the service 
demand is fulfilled. 

In general, this type of schemes requires employing basic 
models for ICI, e.g., lfT4l l22l |2^ . As indicated, another 
approach to CSO is traffic profiling, the case of Gsi m 
M ED, where CSO decisions are taken assuming a cer¬ 
tain knowledge about the service demand ll^ [311 . How¬ 
ever, most of the solutions presented so far employ models 
and assumptions oriented to macrocellular deployments, such 
as III] EH HM . where the use of CSO is not so clear due 
to practical issues ll^ . Other practical aspects, such as the 
number of transitions and handovers are often overlooked 
as well. Recently, CSO has also been studied in the context 
of infrastructure sharing IH and by means of stochastic 
geometry Qa. 

Thus, in the light of these observations, this paper presents a 
novel multi objective frameworl0that 1) includes the strengths 
of previous proposals, 2) overcomes many of their drawbacks, 
and 3) extends the analysis to address aspects that have not 
fully been investigated. 

Main contribution: a framework for multiobjective op¬ 
timization for CSO that explicitly takes into account a 
statistical description of the service demand distribution 
when computing the performance metrics. 

As it will be illustrated by the numerical results, this idea 
has the following advantages: 

/ The use of the spatial service demand distribution (repre¬ 
sented by a spatial probability density function) allows to 

*A preliminary version of this work was presented in fTSl . From it, a US 
patent application has also been made: US patent application no: 14/334,134, 
application date: 17 July 2014. 
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the proposed algorithm to rapidly identify network topolo¬ 
gies, i.e., on/off patterns, providing higher capacity to areas 
where high service demand is more likely to appear, and 
hence, the search space and required computational effort 
is significantly reduced. 

/ Given that, in general, traffic profiles are stable in time 
scales of dozens of minutes, the (computationally-heavy) 
optimization can be done offline and required topologies 
can be applied as needed as these traffic profiles are 
recognized/observed during network operation. This fea¬ 
ture makes the implementation of CSO feasible, given 
that the required BS coordination and real-time processing 
are significantly reduced. However, the conceptual idea 
will still be valid when new paradigms (cloud computing, 
software defined networking, and network virtualization, 
see d) allow faster computation and information exchange 
among network nodes. Thus, more dynamic traffic profile 
recognition and optimization will also be possible under 
the proposed framework. 

/ The proposed optimization formulation allows considering 
several downlink and uplink coverage criteria, such as 
minimum received power and Signal to Interference plus 
Noise Ratio (SINR). 

/ Given that there exist a correlation between the topologies 
that are specific for a given traffic profile, the number of 
handovers and transitions is minimized. 


Finally, the following set of secondary/minor contributions 
are also presented in this work; 

1) Analysis of the impact of the most extended interference 
models (FL and LC) on the performance of CSO. 

2) A quantitative assessment of how CSO operation affects 
the critical uplink power consumption (on user equip¬ 
ment side). To the best of the authors’ knowledge, this 
aspect/criterion has been overlooked in most of previous 
studies, and only recently in M it has been integrated 
within an optimization framework. 

3) While the multiobjective problem formulation presented 
herein can be solved by means of standard stochastic opti¬ 
mization tools ||3^ . an alternative iterative algorithm is also 
proposed for computing the important tradeoff between the 
number of active cells and network capacity. As it will be 
shown shortly, although its performance is slightly inferior 
to stochastic search, it is significantly faster, and hence, it 
could be used when real-time operation becomes feasible. 

The rest of the paper is organized as follows: the next section 
presents the system model. The proposed Multiobjective Opti¬ 
mization (MO) framework (performance metrics and problem 
formulation) is explained in Section m Section describes 
the evaluation setting and benchmarks used in simulations. 
Numerical results are also analyzed therein. Section closes 
the article with conclusions and research directions. 


II. System Model 
A. Definitions and notations 

In this study, an OFDMA cellular network is considered. 
The system bandwidth is B and the network is composed of 
L BSs that can be independently switched off/on. The indices 


of the BSs are contained in the set C = {0,1, • • • ,L—1}. The 
set A, composed of A = |A| small area elements, represents 
the spatial domain to which the network provides service. It 
is assumed that within each area elements a G A, the average 
received power is constant. The maximum transmit power per 
cell is Pmax- The network geometry is captured by the path loss 
matrix G G (distance dependent attenuation, antenna 

gains, and shadowing). The vectors pps and po, both G K^, 
indicate the transmit power at each cell in Pilot Signals (PS), 
and data channels, respectively. Cell selection is based on the 
average PS received power, which can be calculated by means 
of the following expression: 

Rps = G • diag (pps © x), Rps G (1) 

The operator © denotes Hadamard (pointwise) operations. 
The vector x G {0,1}^ indicates which cells are active and 
which ones are switched off. Hereafter, x is also referred to 
as ‘network topology’. Thus, the matrix Rps in Q contains 
the PS received power, i.e., Rps(a, Z) indicates the received 
power from the Z* BS in the a* area element. Of interest is 
the Number of Active Cells (NAC) in each network topology 
as energy consumption is related to it. The a* area element 
(the a* row in Rps) is served by cell Z* if 

Z* = argmax Rps(a, Z). (2) 

iG{0,l,... ,L-1} 

The dependence of Z* on x has not been explicitly indicated for 
the sake of clarity. Based on ([T]| and Q, the binary coverage 
matrices S and S‘^ G can be obtained. If the a* pixel is 

served by Z*, then S(a,l*) = 1. S‘^ is the binary complement 
of S. Hence, the coverage pattern, implicitly defined in S, is 
a function of x. The cell selection rule indicated by Q can 
be regarded as a connectivity function : A —> £ U {—1}. 
If /c(a) = —1, the a* ^ea element is out of coverage, i.e., 

• the received power in a G A (Rps (a,Z*)) is smaller than 

i.e.,Rps(a,Z*)<pRr„, 

• the SINR ('0a) in the area element a G A is smaller than 

i.e., Ipa 0min? t)l 

• the path loss G (a, Z*) between the area element and its 

server is greater than i-®-’ G1 (a,Z*)>GUL^. In prac¬ 
tice, is the maximum path-loss obtained from the 

uplink link budget (a design criterion). 

The cell A; is the subset of A served by the Z* BS. Thus, 
A; = { a G A : /c(a) = Z }, where AiDAj 'i j. The 
set Ac is the subset of A that are associated to one BS. Thus, 
Ac = { a G A : /c(a) —1} = Uig£ The vector T G 

corresponds to the spatial service demand distribution. Thus, 
r(a) indicates the probability, in the event of a new user, that 
the a* pixel has the user on it, and hence, • 1 = 1. It should 
be noted that F is time-dependent, however it is reasonable to 
assume that F is constant during fixed intervals ll30ll . In order 
to represent the service demand volume, two parameters are 
considered; inter-arrival time (A) and session time (p). Both 
are modeled as exponentially distributed random variables. 
Thus, service demand’s spatial distribution and volume are 
described by F and the first order statistics of A and /i, i.e., 
E{A} and E{p,}, respectively. It is assumed that the QoS of a 
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user is satisfied if the target rate (rmin) is fulfilled. Hence, the 
total service demand volume (i?) in A is given by 


where 


ra 


aeA 

• r(o) • rmin [bps] (4) 


corresponds to the average demand in the a* area element. The 
previous model for the service demand can easily be extended 
to the general case of more than one service to account with 
the fact that service time, inter-arrival time, spatial distribution, 
and target rate can be service-specific. Assuming that there are 
Ns service classes, each of them with different characteristics, 
i.e., F'^, and for c = 1,2, ■■■ ,Ns, Q can be 

rewritten as follows: 


r 


s 

a 




[bps]. 


(5) 


The resulting spatial service demand distribution (that is 
required to compute the performance metrics introduced later 
on) can be obtained by considering the resulting demand as 
follows: 

r^(a) = s - (6) 

Hereafter, one single service class (possibly the result of a mix 
of many others) is assumed for the sake of clarity, and hence, 
one single set of parameters (p. A, F, and r^in) are considered. 

Definition II.l (Cell load). The load of the Z* BS (ai(t)), at 
any given time t, is defined as the fraction of the available 
resources (bandwidth) that are being used. 

The average load of the BS is ai = E{ai(f)}. Thus, the 
vector a = [ (Xo Oil • ■ ■ Q^l-i ] indicates the load conditions 
in the network, on average. Note that if x(/) = 0, then 
ol{1) = 0. As the reader can easily infer, as long as ci < 1, 
the network topology (x) is able to satisfy, on average, the 
service demand given by F, E{A}, and E{p}, and hence, it 
can be said that x is adequate. 


Definition II.2 (Network capacity). The network capacity 
(Vcap) is defined as the maximum service demand volume such 
that a < 1. Thus, Vcap — max V : a < 1. 

Definition II.3 (Saturation point). The saturation point (Vsat) 
is the minimum service demand volume such that o; > 1. 
Thus, Vsat — min V : a > 1. 


As indicated, different models can be used for modeling ICI. 
In this work, two models are considered: ‘full load" and ‘load 
coupling’. Recall that in full load, active cells are assumed 
to have full load, i.e., ai — 1, if x(Z) = 1, and ai = 0, if 
x(() = 0. In case of load coupling, the ICI created by each cell 
is proportional to its load. An iterative algorithm to estimate 
the cell load coupling is provided in Appendix Thus, the 
vector fF S representing the average SINR at each area 
element is given by 

= [(S 0 G) ■ (pD 0 x)] 0 [[(S" 0 G) ■ (pD 0 x 0 a)] © cr^] . 

(7) 


The operators 0 and 0 denote Hadamard (pointwise) ope¬ 
rations. By means of average SINR figures as function 
of the network topology (x) are obtained. Since load levels 
also depend on SINR values, the load coupling generates a 
system of non-linear equations which have a unique non¬ 
negative (a > 0) solution 1291 . In order to estimate a, let’s 
look at the average SINR at area element level. The average 
SINR at a & Ai can be expressed as 


V’(a) = - 


Po{l)-G{a,l) 


( 8 ) 


X! • PdO’) • G(a,j) I 0 0-2 

In ([^, the ICI coming from neighbor BSs is proportional to 
their average loads (Aj’s). It is customary to define link perfor¬ 
mance in terms of '!/)(a) by means of a concave (e.g., logarith¬ 
mic) function (/lp) of it, such that = /Lp(^/’(a)) [bps/Hz]. 
The bandwidth requirement of a single user in a & Ai to 
satisfy the QoS can be obtained as 


K{a) = 


[Hz], 


/Lp(V'(a)) 

and the average load {on) in the Z* BS would be given by 

1 


(9) 


Oil = 


B. 


K-bi, 


sys 


where 


and 


N = 


Via) 


E{A}’ 


( 10 ) 


( 11 ) 



K{a) [Hz]. (12) 


In (101, and bi are the average number of users and 
bandwidth consumption in BS I, respectively. 

In order to take into account the coverage criteria and pena¬ 
lize solutions with coverage holes, i.e., a significant number of 
area elements without coverage, the spectral efficiency of the 


a“ area element is stored in the vector H € 


and it is com¬ 


puted according to the following rule: ha = v(a) • /Lp(^/'a)- 
The binary vector v S {0,1}'’^ indicates if the o* is out 
of coverage. Therefore, if the a* area element is in outage, 
v(a) = 1, and 0 otherwise. 

Finally, the list of symbols is provided in Table 


III. Metrics, Problem Formulation, and Solution 
A. Multiobjective optimization: basics 

In order to study the tradeoffs in CSO, the use of mul¬ 
tiobjective optimization has been considered. Multi objective 
optimization is the discipline that focuses on the resolution 
of the problems involving the simultaneous optimization of 
several conflicting objectives, and hence, it is a convenient 
tool to investigate CSO, where the two fundamental metrics, 
energy consumption and network capacity, are in conflict. The 
target is to find a subset of good solutions X* from a set X 
according to a set of criteria iF = {/i,/ 2 ,’’’ !/|j='|}^ with 
cardinality greater than one. In general, the objectives are 
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TABLE II 
Basic notation. 


Symbol 

Description 

Symbol 

Description 

B 

System bandwidth 

Pma.x 

Maximum transmit power per cell 

L 

Number of BSs 

c 

Set with the base stations’ indexes 

A 

Number of area elements 

A 

Set of area elements in the target area 

G 

Path-loss matrix 

Pps, Pd 

Power vectors; pilots and data channels 

R-ps 

Received power matrix 

X 

Network topology 

A 

Coverage area (Ac C A) 

Ai 

Coverage of the Z* BS 

Rps 

Received power matrix 

S, S‘= 

Coverage matrices, i.e., coverage of each cell 

fc 

Connectivity function (cell selection) 

/lp 

Link performance model 

P - 

min 

Minimum received power 

V^min 

Minimum SINR 

^UL 

^max 

Maximum path-loss 

^min 

Minimum target rate (QoS criterion) 

A 

Inter-arrival time 

p 

Session time 

r 

Spatial demand distribution 

OL 

Average load vector 

V 

Service demand volume 


Average SINR vector 

H 

Spectral efficiency vector 

«UL 

Uplink fractional compensation 

V 

Coverage vector 

«cov 

Coverage threshold 

n 

Inverse of cell’s size 




in conflict, and so, improving one of them implies worsening 
another. Consequently, it makes no sense to talk about a single 
global optimum, and hence, the notion of an optimum set X* 
becomes very important. A central concept in multiobjective 
optimization is Pareto efficiency. A solution x* S A’ has 
Pareto efficiency if and only if there does not exist a solution 
X G X, such that x dominates x*. A solution Xi is preferred 
to (dominates) another solution X 2 , (xi X 2 ), if Xi is better 
than X 2 in at least one criterion and not worse than any of 
the remaining ones. The set X* of Pareto efficient solutions is 
called optimal nondominated set and its image is known as the 
Optimal Pareto Front (OPF). In multi objective optimization, 
it is unusual to obtain the OPF due to problem complexity; 
instead, a near-optimal or estimated Pareto front (PF) is found. 
Readers are referred to ||3?I for an in-depth discussion. 

B. Performance metrics 

The following performance metrics have been considerecQ 
• The number of active cells (/i). Under the full-load assump¬ 
tion, energy consumption is proportional to the number of 
active cells ll22l l26ll : 


case, i.e., r(a) = l/A, V a G A. The vector n G 
contains the inverse of the sum of each column in S, i.e., 
the number of pixels served by each cell. It is assumed 
that each user is served by one cell at a time. This vector 
is used to distribute the capacity of each cell evenly over 
its coverage area, i.e., the bandwidth is shared equally by 
the area elements belonging to each cell. This improves the 
fairness in the long run similar to the proportional fairness 
policy that tends to share the resources equally among users 
as time passes. This fairness notion results in decreasing the 
individual rates as the number of users increases. This effect 
is also captured by n as the bandwidth per area element is 
inversely proportional to the size of the cell. 

• Cell edge performance (/a). The 5* percentile of the pixel 
rate Cumulative Distribution Function (CDF) is commonly 
used to provide an indicator for cell edge performance IM). 
A vector with the weighted average rate at area element 
level can be obtained as follows: 


r = A • (H © r) © 


S • (n'^ • diag(S))^ 


(15) 


Then, the percentile 5 is given by 


/i = X • 1. 


(13) 


/3 = r'(0.05-A). (16) 


• Average network capacity (/ 2 ). This metric is based on the 
expected value of the spectral efficiency at area element 
level. Hence, the effect of the spatial service demand dis¬ 
tribution (F) must be considered. The metric is defined as 
follows: 

/2 = (i?-A)[[(H©r)T.S]©n] -1. (14) 

The vector H © F corresponds to the weighted spectral 
efficiency of each area element. The idea is to give more im¬ 
portance to the network topologies (x’s) that provide better 
aggregate capacity {f^) to the areas with higher service 
demand. In ( [T4l l, A (the number of area elements) is used to 
normalize the obtained capacity to the uniform distribution 

^In the definition of some metrics, the dependence with x is not explicit, 
however, it is important to note that all of them depend on x, i.e., the network 
topology. 


The vector r' is a sorted (ascending order) version of r. 

• Uplink power consumption (ff)- In order to provide an 
estimate of the uplink power consumption of any network 
topology, a fractional compensation similar to the Open 
Loop Power Control (OLPC) used in Long Term Evolu¬ 
tion (LTE) is considered ll35ll . It is given by 

• E E r(a) • (Po + f^uL • G (a,()), 

^ ('^1 1=0 aGA 

(17) 

where Pq is a design parameter that depends on the al¬ 
located bandwidth and target Signal-to-Noise Ratio (SNR) 
and KuL G [0,1] is the (network controlled) fractional 
compensation factor. 

• Load dependent power consumption (/o). In order to es¬ 
timate the network power consumption under the load 
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Topology 1 Topology 2 
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Fig. 1. Illustration of network-initiated handover due to CSO operation. 


coupling assumption, the parameterized BS power model 
proposed in ll36l has been used. Thus, 

L-l 

/5 = E^pc(«0, (18) 

1=0 

where ^ function that gives the power consump¬ 
tion in the BS as function of its load. Essentially, in 

this model, there is a hxed power consumption (Pq) that is 
independent of the load but that can be further reduced (till 
Peso) if the base station is switched-off. Moreover, there is a 
part that grows linearly with the load till a maximum power 
consumption (Pmax) that obviously contains the transmitted 
power over the air interface (PjJJax)- 
• Load dispersion (/g). As it will be shown, load dispersion in 
load coupling conditions is an important parameter because 
it measures how well distributed the service demand is. In 
order to quantify this value, the Coefficient of Variation is 
considered. Thus, 


stdjo:} 

meanja} 


(19) 


• Handovers. In the context of CSO, handovers are a quite 
important concern m. Handovers are produced when users 
need to be associated to another base station because their 
serving cells are switched-off. In practice, handovers are 
mainly produced due to users mobility, but independently of 
the type, either user- or network-triggered, handovers require 
a certain time and signaling, both at the air interface and 
core network. Thus, the CSO operation should, as much 
as possible, minimize the number of handover, i.e., the 
transition from one topology to another should be done with 
the minimal impact and/or cost. A pictorial representation 
of the aforementioned situation is shown in Figure [T] Thus, 
handovers are considered herein as an important perfor¬ 
mance metric. 


C. Multiobjective problem formulation 

The multiobjective optimization problem considered herein 
can be formulated as follows: 

optimize f(x) = [fi (x), fj (x)], (20a) 

subject to: 

(A-i-(vT.l)) < «cov, (20b) 

xG {0,1}^, XT^O. (20c) 

Problem ( |20| ) proposes the simultaneous optimization of two 

of the previously introduced performance metrics as follows: 


• Full load: if full load is assumed as model for intercell 
interference, i = 1 and j G {2,3,4} in ( |20al l. 

• Foad coupling: if load coupling is assumed as model for 
intercell interference, i = 5 and j = 6 in ( |20a| i. 

The previous optimization scheme allows to study and 
characterized the tradeoffs between conflicting metrics (see 
Section |III-B| i in a deployment-specific manner. Con¬ 
straints ( |20b| l and ( |20c[ l correspond to the coverage criterion 
and feasible set definition, respectively. 

In general, solving multiobjective problems such as is 
very difficult 1^ . Indeed, ( |20| i is a combinatorial problem 
that belongs to the class NP-complete, and hence, optimal 
solutions cannot be found in polynomial time. The domain 
(search space) dehned by the optimization variable (x, the 
on/off pattern) is a set of size 2^ — 1, where L is the number 
of BSs. The objective space (or image) is defined by the 
objective functions, and due to their mathematical structure, 
it is highly non-linear, non-convex, and full of discontinuities 
and local optima llJTll . Certain algorithms such as Simplex 
are susceptible to be trapped in local optima, while other 
optimization techniques, such as Sequential Quadratic Pro¬ 
gramming 091 . require convexity to guarantee convergence. 
Moreover, traditional constrained optimization, in which only 
one objective function is optimized subject to a set of con¬ 
straints on the remaining ones, has the drawback of limiting 
the visibility of the whole objective space. For this reason, 
heuristic-based algorithms are popular approaches in CSO as 
it was seen in Section but unfortunately, by means of this 
type of solutions it is very difficult to address multi objective 
optimization problems. In order to overcome this difficulty, the 
use of Multiobjective evolutionary algorithms (MOEAs) ll40ll 
is proposed herein as described next. 


D. Multiobjective evolutionary algorithms 

As it was mentioned, heuristic solutions are usually 
problem-specihc and typically used for single-objective op¬ 
timization. Thus, the so-called ‘metaheuristics’ have become 
an active research held Ell. Metaheuristics can be used 
to solve very general kind of multiobjective optimization 
problems, such as the CSO formulation presented herein. 
Indeed, ( |20l i requires a tool able to 1) hnd good solutions 
by efficiently exploring the search space, and 2) operate 
efficiently with multiple criteria and a large number of design 
variables. In addition, it should not have strong requirements, 
such as convexity or continuity. Multiobjective evolutionary 
algorithms (MOEAs) HOl fulhll the previous goals, and hence, 
their use is proposed to deal with the CSO framework pre¬ 
sented herein. MOEAs are population-based metaheuristics 
that simulate the process of natural evolution and they are 
convenient due to their black-box nature that requires no 
assumption on the objective functions. 

Thus, the Nondominated Sorting Genetic Algorithm 
II (NSGA-II) HTl is employed herein to solve ( |20] l. NSGA- 
II is accepted and well-recognized as a reference in the 
held of evolutionary optimization as it has desirable features, 
such as elitism (the ability to preserve good solutions), and 
mechanisms to hexibly improve convergence and distribution. 
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Further details can be found in ia. One key insight for 
selecting evolutionary (genetic) algorithms is that, in CSO, a 
certain correlation is expected among network topologies that 
are suitable for a given spatial service demand distribution, 
i.e., they are expected to be similar, with more cells where 
the traffic is concentrated. The operation in evolutionary 
algorithms precisely does that, i.e., once a good (Pareto 
efficient) solution (network topology) is found, the algorithm 
iteratively try to improve it by 1) combining it with other 
good solutions (crossover mechanism), and 2) adding random 
minor variations to them (mutation mechanism). The complete 
description of NSGA-11 can be found in BTII . As it will 
be shown, the use of MOEAs provides a quite convenient 
approach to CSO. However, depending on the scale of the 
problem, convergence can be slow, especially if computational 
resources are limited. Thus, based on the insight previously 
indicated, and in order to provide additional possibilities. 
Algorithmic is also proposed for solving (finding the set X* of 
Pareto efficient solutions) a particular, yet important, case of 
(20 1 ; when the number of active cells (/i) and the average 
network capacity (/ 2 ) need to be jointly optimized. Given 
that the need for minimizing the number of transitions is very 
important from a practical point of view. Algorithm [C aims at 
finding a collection of network topologies, all with different 
number of active BSs, featuring 1) the minimum distance 
property, and 2) acceptable performance. In this context, the 
word distance refers to the Hamming distance (du), i.e., the 
number of positions in which the corresponding symbols in 
two different solutions are different. In this manner, for two 
solutions Xi and in a set X^^ featuring the minimum 
distance property, c?h (x^,Xj) = 1 ^ |(xi • 1) — (xj • 1)| = 1 
always holds. Initially, Algorithm [C determines the best topo¬ 
logy with 1 active BS (xi) in line 2. Then, in lines 4-14, for 
each successive number of active cells (NAC = 2, ... , L), 
the algorithm sequentially finds the BSs that should be acti¬ 
vated (resulting in the solution Xj), such that 1) the Hamming 
distance with the previous solution Xj_i is one, and 2) the 
function /2 is maximized. Thus, each solution added to 
provides the biggest increment in terms of /2 with respect 
to the one previously added, and only one off/on transition 
is required. It should be noted that, although not explicitly 
indicated. Algorithm indeed optimizes not only the number 
of active base stations and the network capacity, but also the 
number of transitions when moving from one topology to 
another. Thus, more than two objetives are jointly considered. 
The same applies for ( [20l i, i.e., more than two metrics could 
be considered, at expense of an increase in complexity, but it 
should be taken into account that in the context of CSO, a 
Pareto Front in more than two dimensions could complicate 
the implementation. 


E. Conceptual design and implementation 

Figure illustrates the conceptual design of the proposed 
multiobjective framework. The framework relies on having a 
statistical description of the behavior of the service demand 
(in time and space). Thus, by means of different traffic 
distributions (Fa,), the spatial component of the traffic at 


Algorithm 1: Minimum Distance Algorithm (MDA). 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


input : Xi: Ai = {x € A | x ■ 1 = 1}, |Ai| = L. 
output : Amda- a set of L network topologies. 

C* i — 0; Ap4D i — 0; 
xi ^BestBS (Ai) ; 

Xmd “Tmo U {xi}; 
for each j = 2 : L do 
C* ^ 0; 

for each x G Aj | ^^(x, Xj_i) = 1 do 

Cx ^ / 2 (x); 

if Cx > C* then 
C* ^ Cx; 

Xj ^ x; 

end 

end 

Xmd ^ Xmyj U {xj}; 

end 

return A*; 



Fig. 2. Conceptual design of the MO framework for CSO. 


different moments of the day can be captured. These patterns 
can be considered fairly constant during time intervals of small 
duration (tens of minutes or few hours) Il24ll27l . Starting from 
the knowledge of a given T, network analysis and optimization 
based on is done offline. The main idea is that, for 
different demand conditions (spatial distribution and volume), 
different sets of Pareto efficient network topologies can be 
obtained, i.e., for each there is a corresponding A*. 
These sets of near-optimal solutions (A*’s) can be evaluated 
by means of system level simulations (in which several QoS 
criteria, scheduling policies, and ICI models can be considered 
independently) in order to determine which network topologies 
(x G X*) provide the desired level of QoS. Obviously, the net¬ 
work operator may act rather conservatively in this selection 
process as it will be explained in Section IV Moreover, in 
order to allow for semi-distributed implementation, a cluster- 
based operation is encouraged. The benefit of doing so is 
twofold. First, the demand in relatively small areas covered 
by small cells (e.g., pico-cells in a university campus) can be 
characterized easily. Second, the amount of intercell coordina¬ 
tion is reduced compared with the schemes aiming at operating 
in large urban areas. Since demand profiles are stored and in¬ 
dexed at coordinating points in each cluster, the amount of data 
that need to be exchanged (from time to time) is negligible. 
Instead, different clusters (a certain amount of overlapping can 
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(a) Cellular layout. (b) Spatial service demand (F). 




(c) CDF of pixel prob. (r(a)). (d) Power consumption model. 


Fig. 3. Test case scenario. 


be allowed) can also share information in longer time scales, 
so that better decisions can be made in boundary cells. In 
any case, the idea of identifying traffic profiles and applying 
multiobjective-optimized on/off patterns, is compatible with 
novel paradigms that 1) are being considered for 5G {cloud¬ 
networking and virtualization Elia) and 2) would allow for 
more dynamic and centralized operation. In addition, several 
research contributions in the increasingly research held of 
service demand modeling and pattern recognition Il42lj4^ (a 
research problem out of our scope) are appearing, and hence, 
the method proposed herein can extensively be beneht from 
that activity. 


IV. Performance Evaluation 


A. Simulation conditions and parameters 

The simulation setup is based on the assumptions for 
evaluating the IMT-Advanced systems na. The urban micro¬ 
cell (UMi) downlink scenario was chosen. Fig. shows 
the corresponding cellular layout. As it can be seen, the 
network is composed of 37 small cells (radius = 100 m, 
network area « 1 km^). Fig. 3b corresponds to the (irregular) 
spatial service demand distribution (F) used in the numerical 
examples. The Kullback-Feiber distance V with respect to the 
uniform distribution (Tu) can be used as a measure of the non¬ 
uniformity of the spatial service demand distribution as it is 


shown in Figure 3c This setting is perfectly valid to study 


CSO as in this context gains are obtained from the mismatch 
between demand and supply. Indeed, CSO is about hnding the 
smallest network topology that is compatible 1461 enough with 
the service demand to provide the required QoS and coverage. 

In this study, a Single-Input Single-Output (SISO) is con¬ 
sidered. This assumption does not imply any loss of generality 
as long as all cells use the same scheme. The load-dependent 
power consumption model, based on the parameters given 


^ax [dBm] 
18 33 



Cell ID: 0 
(a) Coverage maps. 



^ax [dBm] 

(b) Coverage vs. Tx power. 



(c) Detectable cells. (d) Candidate servers. 

Fig. 4. Coverage aspects: impact of transmit power. 


in 11361 . for pico-BSs assuming an operating bandwidth of 
5 MHz and a maximum transmission power of 30 dBm is 


shown in Fig. 3d 


Dynamic system level simulations are carried based on 
Monte Carlo experiments. The results compile statistics taken 
from 100 independent experiments each of which has a dura¬ 
tion of 5400 s. At each cell, the scheduler assigns each user 
with a bandwidth such that the target rate (rmin) is satished. 
If the percentage of users that obtain a rate equal to rmin is 
greater or equal to the operator-specihc target QoS (Q), then 
the QoS policy is said to be fuffilled. Thus, in order to satisfy 
the maximum number of users, users are sorted based on their 
spectral efficiency and served accordingly. When there is not 
enough bandwidth to satisfy a user, the resource allocation 
ends. The set of parameters used in simulations is provided 
in Table Calibration and complexity aspects of NSGA-II 
are briefly discussed in Section [TV-FI and additional guidelines 
can be found in El mi- The experimentally obtained setting 
is also shown in Table in 


B. Coverage aspects 

The first part of this section is devoted to illustrate some 
coverage aspects and provides insights into the potential 
impact of the transmit power on the performance of CSO. 
Fig. 1^ provides a qualitative perspective. The figure shows 
the size of the maximum coverage (points in which the 
received PS power is greater than Pmin) for the central BSs 
{I — 0) for two different transmit powers (P^ax = 18 dBm 
and Pmax = 33 dBm). For the sake of clarity, shadowing is 
not considered. A quantitative description is shown in Fig. 4b 
which indicates the percentage of the target area (A) that can 
be covered with different values of Pmax- Note for instance 
that, starting from 18 dBm (15% of coverage), P^ax need to 
be increased more than eight times (up to 30 dBm) to double 
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TABLE III 

Evaluation setting and parameters. 


General setting 


Number of cells (L) 

37 (wraparound, omni) 

Carrier freq. 

2.140 GHz 

Bandwidth (B) 

5 MHz 

Max. BS transmit power (-P^^) 

30 dBm 

Pixels’ resolution 

5 X 5 m^ 

Path loss 

M.2135 Umi (45] 

BS’s height 

15 m 

Noise power (cr^) 

-174 dBm/Hz 

Shadowing 

A^(0.4) [dB] 

Cell selection (/c) 

Highest Rx. power 

Rx. power (Pmin) 

-123 dBm 

Small scale fad. 

As in 1451 

Link performance (/lp) 

Shannon’s formula 

Frac. comp, (kul) 

1.00 

Cov. (kcov) 

0.02 

Max. path loss 

163 dB 

Min. rate (rmin) 

400 kbps 

Target QoS (Q) 

97.5% 

User distribution 

According to T 

Traffic model 

Full buffers 

SINR (l/tmin) 

-7.0 dB 

QoS checking interval 

1 s 






Calibration of NSGA-II 




Population size 

100 

Crossover, prob. 

1.00 

Type of var. 

Discrete 

Mutation prob. 

l/L 

Termination crit. 

Hypervolume < 0.001%, ED 




tion III-B (/i, / 2 , /a, and/ 4 ) are provided. Fig. 5a shows the 


Number of active BS (/i) 
(a) Capacity: MOEA vs. MDA. 



Number of active BS (/i) 
(b) Gains: MOEA vs. MDA. 
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(c) Cell edge performance. 

Eig. 5. Multiobjective optimization results. 
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(d) Uplink Tx power. 


has a big influence on the size of the feasible set in ( 201 , i.e., 
the set of x’s for which Constraint |20b| is fulfilled. Hence, the 
impact of Pmax is significant, mainly in low load conditions. 

C. Estimation of network topologies 


First, the results regarding the solution of (20 1 
for the objectives functions introduced in Subsec¬ 


the coverage (up to 30 %), while reaching 60 % of coverage 
requires less than four times the power required for 30 % of 
coverage. Obviously, this depends on the propagation model, 
but the message is that this analysis should be taken into 
account during the design phase of any CSO strategy in order 
to determine appropriate values for Pmax- In the results shown 


in Figs. 4c and 4d all the cells are active and transmit at the 
same Pmax- Fig- 4c indicates the average number of BS that 
can be detected as a function of Pmax (the average is taken 


over the whole coverage area). Fig. 4d shows the percentage 
of the coverage area in which x BSs (servers) are heard with a 
quality (SINR) within X dB below the one of the best server. 
From these results, it becomes clear that the choice of Pmax 


resulting Pareto Front by solving (20 1 , when i = 1 and j = 2 
in ( | 20 a| l, i.e., the joint optimization of the number of active 
BS (/i) and the average network capacity (/ 2 ), by means of 
MOEAs (algorithm NSGA-II) and Algorithmic As expected, 
the use of evolutionary optimization provides better solutions 
than Algorithmic i.e., greater values of /2 for the same value 
of /i. However, it is important to recall that the solutions 
obtained through Alg orithm |C feature the minimum distance 
property (see Section |llI-C| ), and that. Algorithm |C (C>(P^)) 
is, in case of small-to-moderate cluster size, less complex than 
NSGA-II {0{N^ - IPI), N: population size). A quantitative 
perspective of such performance gap is shown in Fig. |^ The 
blue/circle pattern corresponds to the gain in terms of /2 for 
each value of fi indicated in the left vertical axis as ‘Average 
capacity gain’. As a result of the combinatorial nature of 
NSGA-II, the gains are higher when network topologies are 
composed of less BSs, i.e., small values of fi. The red/square 
pattern shows the capacity gain per cell, indicated in the 
right vertical axis. It can be seen that the gain of using 
MOEA is around 1 Mbps/cell in topologies with less than 
20 active BSs (/i < 20). Hence, the use of MOEAs implies 
better network topologies in cases where the computational 
complexity can be afforded. The resulting Pareto Eront by 
solving ( [20| ), for (i = 1, j = 3) and (* = 1, / = 4) in ( |20a| l, 
are shown in Eigs. 5c and |^ respectively. The first case 
illustrates the impact of CSO on cell edge performance. Note 
that while Eig. |^ shows a fairly linear growth of the average 
network capacity with the number of active cells, Eig. |^ 
indicates that cell edge performance (represented by /a) is 
substantially improved only by network topologies featuring 
a higher number of active cells (/i > 27). This result 
clearly suggests that mechanisms for Intercell Interference 
Coordination (ICIC) should be applied together with CSO 
in cases of low load conditions to improve the QoS of cell 


edge users. Eig. 5d illustrates the impact of CSO on the 
power consumption of users (uplink). As it was mentioned, 
the goal is not to determine exact uplink power consumption 
figures, but to create means for comparison among network 
topologies with different number of active BSs. Thus, a 
normalized version of /t (see |T7| is considered. As it can be 
seen, it turns out that the relationship between the number of 
active BSs and the resulting uplink (open-loop-based) power 
consumption is highly nonlinear, being the energy expenditure 
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Fig. 6. Analysis considering cell coupling and load-dependent power consumption. 


considerably high in sparse network topologies (/i < 15). 
Hence, in scenarios where the lifetime of devices should be 
maximized (sensor networks), the use of CSO is not clear. 
Recall that uplink link budget is also considered as a coverage 
criterion. 

To close this subsection. Fig. shows the results corre¬ 
sponding to the solution of ( |20) i for the objective functions 
introduced in Section III-B (/s and /g). According to Defini¬ 
tion 11. 2| and given the spatial demand distribution P (see 
Fig. 3b| , E{A} = 115.0 ms and E{/r} = 119.2 s yield a de¬ 
mand volume {V) equal to Vcap. The resulting load sharing 
patterns (obtained by means of Algorithm]^ for V = Vcap and 
V — 0.5 • Vcap are shown in Fig. Note that increasing V 
results in higher load dispersion. To quantify this. Fig. [&b| 
shows the impact of V on the Coefficient of Variation (CV) of 
the loads (/g). The associated load-dependent power consump¬ 
tion (/g) is also indicated. Note that /g and /g are maximized 
when V = Vsat and V = Vcap, respectively. As expected, the 
load dependent power consumption (/g) is maximized when 
a > 1, i.e., V > Vsat- The dependence of /g on V is 
explained by the strong nonlinearity of ( [T0| ) and the fact that, 
from the load-coupling point of view, ci < 1, and hence, 
no change is expected after V = Vsat- The results shown in 
Figs. 1^ and 1^ are obtained for x = 1, i.e., when all the BSs 
are active. The joint optimization of /g and /g is shown in 
Fig. 6c As it can be seen, there is a conflicting relationship 
between them. The attributes of the extreme solutions (xa 
and xb) in the Pareto Front are indicated. There is also a 
certain correlation between the objectives (/g and /g) and the 
number of active cells (NAC). The topology with the lowest 
energy consumption (/g) requires less active BSs but it has the 
highest load dispersion (/g). Note the difference between the 
highest and lowest loaded BS in xa. In contrast, the best load 
balancing (xb) involves more active BSs, and hence, worst 
values of /g. A comparison among solutions obtained through 
each ICI model, FL and LC, is provided next. 


D. System level simulations 

As indicated earlier, solving ( p0| ) results in a set of Pareto 
efficient (nondominated) network topologies that are specific 
for either a spatial service demand distribution (Afl: full¬ 
load) or a service demand conditions, i.e., spatial demand 


distribution plus volume (Alc^ load-coupling). Recall that 
Afl is obtained by joint optimizing /i and /2 in (20i for a 
given spatial demand distribution (F), while obtaining A’lc 
involves the joint optimization of /g and /g in (j2^ for a 
given r and V (volume). Note that, the ‘full-load’ analysis 
is volume-independent, and hence, it does not require spec¬ 
ify V (full load is assumed for the active cells). Thus, in 
order to evaluate these solutions by means of system level 
simulations, it is initially assumed that at each QoS checking 
interval (evaluation parameters are shown in Table [nl| , the 
(nondominated) network topologies of each set (Afl and Alc) 
are all applied and evaluated. The goal is to create QoS 
statistics for each network topology and load condition. Then, 
the network topology that is able to provide the desired QoS 
(Q % of users are satisfied Q% of time) is selected and 
applied (as indicated in Subsection III-E| i. The comparative 
assessment is shown in Fig. |7] where the legends indicate the 
set the applied network topology belongs to (Afl or Alc) and 
the ICI model (FL or LC) used in the system level trials. 
Fig. 1^ shows the load-dependent power consumption of each 
network topology. Clearly, from the CSO point of view, the 
topologies in Alc result in lower power consumption as they 
feature less active BSs (NAC is indicated in green boxes) given 
that the load-coupling model predicts better SINR than full¬ 
load (see Fig. |^, and hence, network capacity is favored. 
However, as V increases, both models become somehow 
equivalent as the loads tend to 1; as a result, the energy 
consumption is quite similar. Figs. and show the QoS 
level (in terms of the number of satisfied users) that is obtained 
with the selected solution of each set for V = 0.2 • Vcap 
and V — 0.6 • Vcap, respectively. The results make evident 
that the performance of the network topologies in Alc is 
severely degraded if the ICI levels become higher than the ones 
from which they were calculated for, see Alc@FL (full ICI). 
Indeed, the performance of these solutions is sensitive to vari¬ 
ations from the mean values (that happens when considering 
snapshots) in moderate-to-high load conditions, even when the 
load-coupling based ICI is considered, as seen in Fig. 7c for 
Alc@LC. On the other hand, the network topologies in Afl 
provide consistent performance when they are evaluated under 
full load (Afl@FL), and obviously, provide an even better 
performance under load-coupling (Afl@LC) for both demand 
volumes. Hence, given that the energy consumption gain is in 
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Fig. 7. Comparative analysis: Load-coupling vs. Full ICIC MO optimization. 



(c) y = 0.6 • ycap. 





Fig. 8. Performance comparisons. 


(c) Key performance indicators. 


the order of 10 % in the best case, it can be concluded that 
the full load model provides a competitive and somehow safer 
energy-saving vs. QoS tradeoff in the context of CSO. The 
proposed CSO scheme can use either approach. Summarizing; 

/ The MO for FL, i.e., fi and /2 in ( |20l l, is volume- 
independent; offline system level simulations are required 
for each load condition (V), and energy saving is smaller 
in comparison to LC. 

/ In MO for LC, i.e., /s and /g in ( [20| ), is volume 
dependent; different offline optimization procedures are 
required for each load condition (V), and energy saving 
is larger in comparison to FC. 

E. Performance comparison 

In order to provide a wide perspective of the merit of the 
CSO framework presented herein, several recent/representative 


CSO schemes have been used as baselines. Obviously, an 
exhaustive comparison is not feasible. However, the idea is 
to illustrate some pros and cons of different approaches 
and the impact of some design assumptions. The following 
benchmarks are considered; 

/ CeJJ zooming: It was proposed in ll26l . The idea is to 
sequentially switch-off BS starting from the lowest loaded. 
The algorithm ends when a cell cannot be switched-off 
because at least one user cannot be re-allocated. 

/ Improved cell zooming: This scheme is presented in 
and it is similar to the one in f26\ . but it includes a more 
flexible termination criterion that allows to check more cells 
before terminating, and so, more energy-efficient topologies 
can be found. 

/ Load-and-interference aware CSO: The design of this CSO 
scheme presented in El takes into account both the re¬ 
ceived interference and load of each cell to create a ranking 



















































































12 


that is used to sequentially switch-off the cells whose load 
is below a certain threshold. 

/ Set cover based CSO: The CSO scheme proposed in iflTll 
relies on the idea of switch-on BS sequentially according to 
a certain sorting criterion. In this work, the sorting criterion 
is based on the number of users a cell can served in the 
SNR regime. 


The performance comparison is shown in Fig. To make the 
comparison fair, the full-load ICI conditions are considered. 
Figs. 8a and 8b show the average number of active cells and 
QoS (for different service demand volumes), respectively. As 
it can be seen, the best energy saving is obtained by ini, 
although at the expense of QoS degradations. This is due to 
the fact that in ini, users can be easily put in outage. In 
contrast, CSO schemes such as iflTll22ll26l provide the desired 
QoS (as long as C < Vcap) since CSO decisions require 
associating all users. However, this results in an increment 
in the average number of active cells with respect to ini. 
The schemes labeled as ‘MDA’ and ‘MOEA’ correspond to 
the (infeasible) dynamic selection of network topologies from 
the sets Xpl obtained through Algorithms [T] and NSGA-II, 
respectively, which are shown as reference. The performance 
of the proposed MO CSO is indicated by red boxes and 
labeled as ‘MO CSO’. As it can be seen, the proposed scheme 
provides an excellent tradeoff between the required number of 
active cells and the obtained QoS, especially when V > Vcap 
where the performance (QoS) of other CSO is compromised. 
However, the most significant enhancement in the proposed 
scheme is its feasibility. Fig. shows four performance 
indicators: transitions, handovers, QoS, and NAC. Given that 
the network topologies are calculated offline, they can be 
evaluated extensively by means of system level simulations 
(under a wide range of coverage criteria and conditions) to 
further guarantee their real-time performance, i.e., the operator 
can select topologies with more active cells rather than the 
ones which strictly need to guarantee QoS. Therefore, the 
selected network topologies can be applied (without real-time 
complexity) during periods of time in which service demand 
is described by F; as a result, no transitions or handovers 
are induced due to CSO. Hence, feasible yet effective CSO 
performance is achieved. As it was shown earlier, the proposed 
framework is generic, flexible, and no assumption are made 
in regards to, for instance, the cellular layout or objective 
functions; as a result, the framework is also suitable for small¬ 
cell deployments where irregular topologies and heterogeneous 
demand conditions are expected. 


E Complexity analysis and calibration aspects 

To close this section, a complexity overview of the optimiza¬ 
tion algorithms is provided. According to El, the complexity 
of NSGA-II is 0{M ■ N^), where M and N correspond to 
the population size and the number of objective functions, 
respectively. In our case, N = 2 and M can be set depending 
on the scale of the problem. However, there is a consensus 
about the size of the population when using genetic algorithms, 
such as NSGA-II, and it is considered that during calibration 
populations of 20 up to 100 individuals can be used. Values 


greater than 100 hardly achieve significant gains and the same 
global convergence is obtained HTl . Regarding Algorithm 
it’s complexity is 0{LP'), where L is the number of cells 
in the network. In practice, Lp' M ■ N'^ which is a 
significant reduction in terms of complexity that comes at 
expense of some performance. In evolutionary algorithms, a 
termination criterion is usually defined/need. One metric used 
to measure the level of convergence is the the hypervolume 
indicator HOll . It reflects the size of volume dominated by the 
estimated Pareto Front. In this work, the seacrh is terminated 
if the improvement in the hypervolume is smaller than a 
threshold (0.001%) after a certain number of generations (in 
this study, 20). Finally, crossover and mutation probabilities 
are set to 1 and 1 jL (one mutation per solution, on average). 


respectively, as indicated in Table III 


V. Conclusions and Research Directions 

CSO is a promising strategy that allows significant energy 
saving in cellular networks where both radio access net¬ 
work (capacity supply) and service demand are heterogeneous. 
In this article 1) CSO has been carefully analyzed considering 
coverage criteria, ICI models, and practical aspects, such 
as network-initiated handovers and on-off/off-on transitions, 
and 2) a novel MO-based CSO scheme has been introduced. 
The proposed solution succeeds in minimizing the number of 
transitions and handovers caused by the CSO operation and 
it is able to operate without need for heavy computational 
burden as the core processing is done offline. In addition, a 
cluster based-operation have been proposed to allow for semi- 
distributed implementation. The results show that, when com¬ 
pared with previous proposals, the proposed solution provides 
competitive performance in terms of QoS and energy saving 
while offering clear advantages from the feasibility perspective 
as it reduces the number of handovers and transitions. The 
results also highlight the importance of considering coverage 
criteria (in downlink and uplink) and pay attention to the 
selection of operational parameters, e.g., the power allocated 
to PS (typically used as criterion for coverage). 

A comparative analysis between ICI models (full-load and 
load-coupling) indicates that the full-load assumption is a 
safe approach in the context of CSO as it provides natural 
protection against deviations from average load values that are 
1 ) used as input of the algorithm, and 2) inherent of real time 
operation, i.e., discrete realizations of users. The impact of 
CSO on the power consumption of UE has also been studied. 
The results indicate that sparse topologies (few active BSs) 
have a significant impact on uplink power consumption, and 
hence, CSO is not suitable for scenarios with energy-sensitive 
devices such as sensor networks. 

Research on topology adaptation has still a long way until 
its maturity. Eeasible and effective techniques for traffic pattern 
recognition to complement CSO are still in infancy. It is our 
strong belief that CSO, as a promising approach to greener 
networks, is a key piece of a more general set of capabilities 
that will appear in 5G networks, also including promizing 
and disruptive concepts, such as Downlink Uplink Decou¬ 
pling (DUDe). DUDe, where user equipment can transmit 
and receive to and from different base stations, is indeed, a 
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Algorithm 2; Iterative approximation of cells load. 

Inputs : All relevant information and e (termination). 

Output : a £ Load vector. 

1 a® ■<— 1; e ■<— e • 1; fc ■<— 0; /* initializing */ 

2 repeat 

3 fc t— fe + 1; 

4 a* ■<— 

5 for 1 = 1 : |£| do 

6 af 3-Load ( : (Vj G £) A (j 7 ^ ()}) ; 

7 3—df; /* Fast update */ 

8 end 

9 3 — [dg • • • dj,_i]; /* Update: iteration k */ 

10 until e > 

11 return a*; /* Return estimated load vector */ 


clear research direction from the perspective of CSO, where 
both uplink and downlink could be considered as independent 
networks. 


Acknowledgment 


The authors would like to thank Tamer Beitelmal from 
Carleton University, and Dr. Ngoc Dao from Huawei 
Canada Research Centre, for their valuable feedback. This 
work was supported by the Academy of Finland (grants 
287249 and 284811). Mario Garcia-Lozano is funded by 
the Spanish National Science Council through the project 
TEC2014-60258-C2-2-R. 

Appendix A 

Iterative approximation of cell loads 

In order to estimate the average load vector (a). Algorithm]^ 
is proposed. Basically, the estimation of the average load at 
each cell (d;) is refined through each iteration comprising 
Lines[^to|^ In line|^ the function Load () estimates each d, 
based on ((TO]), from the values of previous iterations (where 
d/ = minjl, d/}) and the ones that have been just updated in 
the current iteration (this/ait update is done in line0. Fig.jg 
illustrates the motivation and performance of Algorithm 1^ 
Basically, the use of load coupling provides a more accurate 
estimation of ICI levels in the network as shown in Fig. 
Note that the use of Full ICIC represents a more conservative 
approach. Fig. 9b shows that Algorithm only requires few 
iterations to converge and that this depends on the starting 
point, but in any case convergence is fast. 
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